extern integer isamax_(integer *n, real *sx, integer *incx);
